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Abstract 

We compute moments of non-diagonal correlators with two massive quarks. Re- 
sults are obtained where no restriction on the ratio of the masses is assumed. Both 
analytical and numerical methods are applied in order to evaluate the two-scale 
master integrals at three loops. We provide explicit results for the latter which are 
useful for other calculations. As a by-product we obtain results for the electroweak 
p parameter up to three loops which can be applied to a fourth generation of quarks 
with arbitrary masses. 

PACS numbers: 12.38.Bx, 14.65.Dw, 14.65.Fy 



1 QCD corrections to current correlators 



In recent years many multi-loop results to current correlators became available. In the 
massless approximation the state-of-the art are four-loop integrals which have been applied 
to t and Z-boson decays in order to extract a s [1,2]. On the other hand, in the limit of 
small external momentum the so-called moments have been obtained also up to four-loop 
order [3-11] and both a s and the charm and bottom quark masses could be extracted by 
comparing to experimental data or lattice simulations [12-16]. The currents involved in 
these calculations only couple to one quark flavour. 

As far as the non-diagonal correlators, where the external currents couple to two quark 
flavours with different masses, are concerned results for the moments up to three loops 
are known for the limit where one mass is zero [17-22]. At four- loop order the vector and 
axial-vector correlators have only been computed for vanishing external momentum (see 
Ref. [23,24]) in order to obtain four-loop corrections to the electroweak p parameter in 
the Standard Model [23-25]. 

It has been argued in Ref. [26] that it would be desirable to have at hand moments of the 
heavy-light currents also for general values of the two quark masses. The reason to not 
only consider physical mass values of the involved quarks is related to lattice simulations 
where often a variety of different masses are considered in order to be able to extrapolate 
to the mass values of physical interest. In Ref. [22] this task has been solved with the help 
of expansions around the equal-mass case and around the limit where one mass is much 
smaller than the other. In this paper we present exact results up to three loops and thus 
improve the findings of [22]. 

We adopt the notation from Ref. [22]. Let us nevertheless present the relevant formulas 
in order to have a self-contained presentation of the results. 

The momentum-space correlator formed by a vector (v), axial- vector (a), scalar (s) or 
pseudo-scalar (p) current is given by 

HV + ») nV) + »ni( g 2 ) = iJaxe**(o\Tjfe)$(o)\o), 

q 2 Yi\q 2 ) = z|dxe^(0|T/(x)/t( )|0), (1) 

with 

3l = ^i7/^2 , 3l = ^i7 M 7 5 ^2 , f = ^ifa , f = V^7 5 ^2 • (2) 

In our calculations we are allowed to use anti-commuting 75 since for ip\ ^ ip2 only non- 
singlet diagrams contribute. 

We work in rij-flavour QCD with two massive and nf — 2 massless quarks. In our final 
results we adopt the MS renormalization scheme both for the parameters (strong coupling 
constant a s and the quark masses mi and m 2 ) and the overall renormalization of the 
correlator which in the following is indicated by a bar. 
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In the limit q 2 — > the quantity U s (q 2 ) (and analogously U s L (q 2 )) can be cast in the form 



16vr 2 

n>-l 



where the dimensionless variables 



i 



^ A 7712 (A\ 

z = — =■ , and x = — , (4) 
mf mi 

have been introduced. Note that we assume < x < 1. Results for rri2 > m\ are easily 
obtained by interchanging m\ and m 2 . We refer to the coefficients C^(x) also as moments. 
Their perturbative expansion can be written as 

^ = ^ + -^+(-) 2 ^ + ---, (5) 

7T V 7T / 

where the arguments x and /i (renormalization scale) are suppressed. Note that we have 
the symmetry relations C^(x) = C%(—x) and C%(x) = C®(—x) (see, e.g., Ref. [27]) which 
we use to cross check our results. A further check is provided by computing the moments 
of the longitudinal contribution Il^' a (g 2 ) and compare the result of the moments of U s,p (q 2 ) 
with the help of the relations 

Ctn = (1 + ^+1- (6) 



For completeness, let us mention that the QED-like normalization of the polarization 
function with II 5 (0) = can be obtained from IT 5 (g 2 ) with the help of 

nV) = nV) - ^ (<?o 4 + c -f) . (7) 

Sample Feynman diagrams occurring at one-, two- and three- loop order are shown in 
Fig. 1. 



2 Master integrals for non-diagonal correlators 

For the evaluation of the Feynman diagrams a well-tested chain of programs has been used 
which works hand-in-hand in order to avoid error-prone manual interactions. Since the 
reduction to master integrals is nowadays pretty standard we restrict ourselves to a brief 
description up to this step. The Feynman diagrams are generated with QGRAF [28] and 
the various diagram topologies are identified and transformed to FORM [29, 30] notation 

1 In the following we do not write rrii f \[i) and m£ ' (fi) but suppress the dependence on the renor- 
malization scale \i and the number of flavours rif. 
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Figure 1: Sample diagrams contributing to IT* 5 at one, two and three loops. The thick 
(upper) and thin (lower) lines correspond to quarks with mass m\ and m 2 , respectively, 
and the curly lines represent gluons. 

with the help of q2e and exp [31,32]. In a next step appropriate projectors are applied, 
the expansion in the external momentum is performed and traces are taken. Afterwards 
the scalar integrals are mapped to functions which are then passed to FIRE [33] where the 
reduction to master integrals is performed. 

The master integrals which are needed for our calculation are shown in Figs. 2 and 3. At 
one- and two-loop order there are one and three master integrals, respectively, while at 
three-loop order the reduction leads to three single-scale (c.f. Fig. 2) and 16 two-scale 
master (c.f. Fig. 3) integrals. All one- and two-loop integrals and the single-scale three- 
loop integrals in Fig. 2 are available in the literature (see, e.g., Ref. [34,35]). For the 
three- loop integrals in Fig. 3 only partial results exist. In this work these results are 
checked by independent calculations and results for all integrals are presented. Since the 
problem at hand is symmetric under the exchange of mi and most master integrals 
come in two variants which are marked by the subscripts "a" and "b" . 

2.1 Differential equation method 

Most of the master integrals can be calculated with the method of differential equations 
(see, e.g., Refs. [36-41]). In the following we exemplify the method for the two- loop master 

(2) 

integral Jf of Fig. 3 which is given by 

where d = 4 — 2e is the space-time dimension. The differential equations are derived 
by taking the derivative with respect to x and reducing the resulting integrals to master 
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Figure 2: One-scale master integrals where solid and dotted lines denote massive and 
massless propagators, respectively. 



integrals. This leads to 



X 



dx 1 v ' 1 — x 2 1 v ' x(l — x 2 ) 

where K^ 1 ' (x) denotes the one-loop integral in Fig. 2 with mass xui\ which is given by 



(9) 



K^\x)=M 



(d-4)(d-2) 



(mi)' 



with the normalization 



J\f 



^] tir d / 2 T(3-d/2). 



(10) 



;ii) 



To solve the differential equation it is necessary to know the initial condition at either 

(2) ('9 s ) 

x = or x = 1, which corresponds to the single-scale integrals K{ and K\' shown in 
Fig. 2. In this simple example it is possible to solve the differential equation without 
expanding in e with the result 



jf)( x ) =M 2 m\ 



2(1 -x 



2M-2e 



r( 



2e] 



r i 



x 



2-2e, 



1 — X 



2M-2e 



-l + 6 2 e 



2 C 2 



2 F 1 (2-2e,l-e,2-e,a; 2 ) 



(12) 



This result agrees with the one given in Ref. [34]. In general it is not possible to obtain 
a closed solution as in Eq. (12). In those cases one has to expand the original differen- 
tial equation in e and thus obtains simpler differential equations for the corresponding 
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Figure 3: Two-scale master integrals where thick, thin and dotted lines denote heavy, light 
and massless particles. A dot on a thick line indicates that the corresponding propagator 
is squared and a cross marks propagators which are raised to power minus one. 
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coefficients. In our example one gets from Eq. (9) for the first three coefficients 



2x 



.r- 



2x 



x 2 — 1 



2x 



x- 



4%(x) - 



x z 



■ -AO 



2x(l - 2 log x) 



1 



x- 



2x 



■J 



x- 



(2), 
1,0 



4x t(2 ) , , 2x(l - 21ogx + 21og 2 x) 

-Jl'Ax) + 



(13) 



where the j[ 2 l(x) denotes the coefficient of the n th order in the expansion in e. One 
observes that the structure of the homogeneous part remains the same whereas the in- 
homogeneous contribution becomes successively more complicated. While solving the 
differential equations it is advantageous to use Harmonic Polylogarithms (HPLs) [42], 
which form a special class of functions suitable for these kind of problems. It is par- 
ticularly convenient to use the implementation in Mathematica from Ref. [43,44] which 
allows the application of the method of the variation of constants in order to solve the 
differential equation in Eq. (13). For our example 2 at hand the result looks as follows 



J {2 \x)=M' 



2 m 2 



> 2 + 1 )) + 7 



(1 + x 2 ) + 2x 2 H ( 



x 



+ 



7T 2 Qr 2 -l)-21(l + x 2 ) 
6 



+ 6x 2 H (x) - 4x 2 H , (a;) 



- 2 (1 - x 2 ) (H 1>0 (x) - H_ 1>0 (x)) + 0(e) 



(14) 



which can also be obtained by expanding Eq. (12) and agrees with Ref. [34]. 



Some of the three-loop master integrals with two masses ( J§ a , J^ b , Jg a , J^ b ) are products 
of one- and two-loop integrals and can therefore easily be calculated with the result 

4 3 J(x)=[K^(l)] 2 K^(x), 
4!\x) = kV(1)[ K V(x)} 2 , 

j£\x) = 4 2 \x)kW(x), 

where K^\x) and j[ 2 \x) are given in Eqs. (10) and (14), respectively. 

2 Eq. (13) can be expressed in terms of HPLs with the help of H (x) = log 2 and Ho.o(x) = (log.T) 2 /2. 
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O o t(3) t(3) r(3) ,(3) r(3) j .(3) 

The method of differential equations described above can immediately be applied to the 
the three-loop integrals j[ 3 \ J 2 3 \ jf \ jf \ and J®. j[ 3 ^ and J 2 obey a system of 
coupled differential equations, as do Jg and j[ 3 \ The integrals and jf^ which are 
symmetric under interchange of mi and 777.2 can be obtained by solving a single differential 

equation. Solving the differential equations leads to analytical results in terms of HPLs. 

(3) 

The integrals needed as initial conditions are shown in Fig. 2. The result for J] is given 
in Ref. [41]. In this reference also the result for an integral related to J 2 by integration- 
by-parts identities is listed. We find agreement including terms of order e. 

To our knowledge the results for J% 3 \ jf\ and J® are new. They read 



A 3) (x) =Af 3 mf 



x 

3^ 



+ h (12 ( _1 + 16x2 ~ x ^ ~ x2H °^) 

+ i (~ (3 - 16x 2 + 3x 4 ) + l -x 2 (-8 + x 2 ) H (x) + 2x 2 H , (x)^ 

+ — (-145 (1 + x 4 ) + 4vr 2 (-1 + x 4 ) + 8x 2 (35 + 4Cs)) 
48 

+ -^x 2 (-4 (30 + 7T 2 ) + 45x 2 ) H (x) - x 2 (-8 + x 2 ) H , (x) - 4x 2 H 0i0 , (x) 



- (1 - x 4 ) (H li0 (x) - H_ 1)0 (x)) - 4x 2 (H 2 , (x) - H_ 2 , (x)) + 0(e) 



(16) 



jf(x)=Af 3 ™ 6 



777; 



3e J \36 v 7 



x 



- ( — (-10 + 69x 2 - 26x 4 + 5x 6 ) 
c \ 24 



+ x 2 ^-4 + x 2 - j J H (x) + 2x 2 H ,o(x) 
+ 1 (-145 - 4tt 2 ) + 1 (203 - 4tt 2 + 32C 3 ) x 2 



37 7T 2 

T + 6 



X 



+ — (145 - 4vr 2 ) x 6 - ^x 2 (124 + 4tt 2 - 82x 2 + 15x 4 ) H (x) 

+ ^x 2 (24 - 6x 2 + x 4 ) Ho l0 (x) - 4x 2 H , ,o(x) - 4x 2 (H 2 , (x) - H_ 2 , (x)) 



- - (2 + 3x 2 - 6x 4 + x 6 ) (Hi, (x) - H_ 1)0 (x)) + 0(e) 



(17) 



J?J(x) =Af 3 ml 



-1 - x' 



^(-2-^ + 2x 2 H (x; 
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+ (-25 - 17x 2 + 7T 2 (-1 + x 2 )) + 10x 2 H (x) 
-4x 2 H 0i0 (x) - 4 (1 - x 2 ) (Hi, (x) - H_ 1)0 (x)) 

+ - (-90 + 5tt 2 (-1 + x 2 ) + 22C 3 - 7x 2 (7 + 2C 3 )) + 34x 2 H (x) 

tt 2 Cl — 4r 2 4- 3r 4 l 

- 20x 2 H ,o(x) + 8x 2 H 0A0 (x) + — ^ q 2 ( H i( x ) ~ H -i( x )) 

- 20 (1 - x 2 ) (H 1)0 (x) - H_ 1)0 (x)) + 8(1- x 2 ) (H 1A0 (x) - H_ 1A0 (x)) 
+ 4 f-4 + \ + 3x 2 ) (H 1)1)0 (x) + H_! _ 1)0 (x)) 



x 2 



4 



16 - — - 12x 2 ) (H li _ 1:0 (x) + H_ lilj0 (z)) + 0(e) 



J) b can be obtained from Jj a by interchanging mi and m%. We have performed explicit 

(3) (3) 

calculations for J^ a and J) b and have used the symmetry relation as a check. Note that 
the HPLs exhibit cuts along the positive real axis and thus one has to be careful to use the 
proper analytic continuation. The analytic results for all master integrals can be found 
in the file TwoMassTadpoles .m which can be obtained from Ref. [45]. 

There are more checks to verify the obtained results. First of all we have checked that the 
results for the master integrals satisfy the original differential equations. Furthermore, 
since it is sufficient to use the value of the integral at x = or x = 1 as initial condition 
for the differential equation the value at the other boundary can be used as cross check. 
A further check constitutes the successful comparison to the expansions around x = 
and x = 1 [22]. 



4 t(3) t(3) .(3) , .(3) 
J 8a , j gb , j Qa ana j % 

The integrals jj^ and jj^ (and correspondingly jj^ and Jgf) obey a system of two 
coupled differential equations which could not be solved analytically using the method 
described above. Providing initial conditions a numerical solution is possible, however, 
the achieved accuracy for the master integrals is not sufficient to compute higher order 
moments of the current correlators since there are large numerical cancellations between 
contributions from different master integrals. Thus we decided to apply the Mellin-Barnes 
method [46,47] which provides for a given value of x a high-precision numerical result with 
about 30 significant digits. Considering the integral and allowing for generic indices 
n\ and n 2 on a line with mass mi and m 2 , respectively, one can derive the two-dimensional 
Mellin-Barnes representation 

dz! / dz 2 m- 2(ni+n2+22+3e - 4) m 2z2 7r 4 - 3 T {- Zl ) T (-z 2 ) 

-ioo J —ioo 
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r (-21 - e + 1) r (-m - z x - 2e + 3) r (-m - *i - e + 2) 

r (— n 2 - z 2 - e + 2) T (m + n 2 + Z\ + z 2 + 2e - 3) 
T (m + n 2 + ^ + z 2 + 3e - 4) / [4I\2 - e)r (m) T (n 2 ) 
r(-m-2^-2e + 3)] , (19) 

where the integration contour is chosen to separate the poles from the Gamma functions 
of the form T(z\ + z 2 + . . .) from the ones originating from r(— Z\ + . . .) or r(— z 2 + ...). 
We used the package MB [48] to evaluate the Mellin-Barnes representation numerically for 
0.1 < x < 1 in steps of 0.005 with high precision. 3 The results can be interpolated to 
obtain results for all values of x. An important cross check of the numerical approach 
constitutes the comparison of the exact analytical results for the poles which can be 
obtained from asymptotic expansions or the Mellin-Barnes representation. They are given 
by 



(3) 



Af 3 m 



1 



(1 + x 2 ) + 



1 



15 

T 



Ax 2 - 



x 
12 



3x 2 log(x' 



1 /65 
f - ( — + 10x 2 



-x 



-12x' 



r 



log(x) + 3x 2 log 2 



x, 



■I 



(3) 



C(e c 



AA 3 mf 



_L + 4x 2 + ^ 



— 3x 2 log(x) 



1 /— 5 f,^ 45 

-6x 4 log(z)) + - — + 10x 2 + —x 4 - (12x 2 + —x 4 ) log(x) 
e \ 8 8 2 



+ (3x 2 + 18x 4 )log 2 (x)) +C(e c 



J, 



(3) 



9d 



AA 3 m? 



2 _ x 2 

3 ~ ~3 



1 



3 
2 



5x 
~6~ 



(3) 



i)b 



Af 3 m 2 



5 4x 2 
3 ~ ~ 

.2 j 



— X 



+ 2x 2 log(x) — x 2 log 2 (x 
5x 2 



+ x 2 log(x^ 



O(e ) 



6 



3x log(x y 



+ - (1 - 4x 2 + 6x 2 log(x) - 3x 2 log 2 (x)) + O(e ) 



Note that for these integrals the order e 1 terms are not needed for the final results. 



(20) 



3 For x < 0.1 we can safely rely on expansions, see also Section 2.5. Actually, for the practical 
luation of the master integi 
the one around 1 for x > 0.5. 



(3) (3) (3) (3) 

evaluation of the master integrals Jg a , Jg b , Jg a and Jg b we use the expansion around for x < 0.2 and 
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2.5 jgl and J® 

(3) (3) 

The differential equations for J 10 ^ and J{ Qb only contain inhomogeneous parts involving the 
integrals jj^ , and , J$l\ respectively, and can therefore not be solved analytically 
However, in this case the straightforward numerical solution of the differential equations 
would lead to results which are sufficiently precise for the moments considered in this 
paper. 

A more precise result for J^ a and j[^ b which is furthermore simpler to handle can be 
obtained with the help of expansions around x = and x — 1. In contrast to Ref. [22] 
where expansions of the whole correlator have been considered and thus the expansion 
depth was limited to 8 and 9 terms, respectively, in this paper results up to order x 38 and 
(1 — a;) 28 could be obtained for the scalar integrals J^ a and J^. 4 The upper plots in 
Fig. 4 show the results for J$ a and J^ b using the expansions around x = (solid) and 
x = 1 (dotted). The difference of these results is shown in the bottom plots. One observes 

(3) (3) 

perfect agreement over a wide range of x for Ji 0a . In the case of J^ ob good agreement is 
only found for 0.2 < x < 0.4 which suggests to use the x — > results for x < 0.3 and 
the x — > 1 expressions for x > 0.3. This has been implemented in the program for the 
numerical evaluation of the master integrals. 

(3) (3) 

The divergent parts of Ji 0a and Jl ob can be obtained in analytical form with the help of 
an asymptotic expansion. They read 

4H = A/- 3 ^ + 0( e °), (21) 
J® = Af 3 ^ + O(e ). (22) 

Also for these integrals the order e 1 terms do not enter the final results. Thus numeri- 
cal expressions are only needed for the finite contributions to the moments whereas the 
cancellation of the poles can be checked analytically. 

All available information about the master integrals is provided in the above mentioned 
Mathematica file TwoMassTadpoles .m [45]. 



3 Results 
3.1 Moments 

In this Section we present numerical results for the first four moments of the vector, axial- 
vector, scalar and pseudo-scalar current correlators. We adapt the colour factors to 5*^7(3) 
and set the number of massless quarks to three. Furthermore, we set the renormalization 

4 Note that there are cancellations when constructing the physical moments using the expansion around 
x = 1. Thus some of the moments are only known up to order (1 — x) 15 . 
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Figure 4: Comparison of asymptotic expansions for j[^ a (left) and j[^ b (right). The plots 
on the top show the O(e ) results for the master integral where the dotted and solid curve 

(3) 

corresponds to the expansion around x = 1 and x = 0, respectively. In the case of J{ Qa 
the dotted curve is barely visible since it is on top of the solid one. The plots on the 
bottom show the difference of the two expansions. 



scale to rrii. More general results that still contain the invariants of SU(N C ) and labels 
for the light and heavy quarks can be found in the Mathematica file coefhl2.m [45]. 

For comparison we also show in Figs. 5 to 8 the results obtained in Ref. [22] as dashed 
lines. In Ref. [22] expansions around x = and x = 1 have been computed which showed 
good convergence properties up to x 0.1 and x ;> 0.5, respectively. Inbetween an 
interpolation has been performed. 



In most cases good agreement is found. Small deviations for x ~ 0.2 



for Ci 2 ^' v and, to a lesser extent, also for C 



(2),v £f(2),a ^(2),a 



and C, 



. . 0.3 are observed 
They demonstrate 



the importance of the calculation performed in this paper. 

From the webpage [45] it is possible to download a Mathematica file which provides the 
results for all moments of Figs. 5 to 8. Furthermore, also the results for Cn with n — — 1 
and n = and for CI n and C*£ n with n = — 1, . . . , 4 are included. The one- and two- loop 
expressions are available up to n = 9. 



12 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



x x 

Figure 5: Dependence on x of the first four moments of the vector current correlator. The 
dotted lines correspond to the interpolation results obtained in Ref. [22]. 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



x x 

Figure 7: Dependence on x of the first four moments of the scalar current correlator. The 
dotted lines correspond to the interpolation results obtained in Ref. [22]. 
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3.2 p parameter 



In this section we discuss an immediate application of G 5 _ x of the vector and axial-vector 
current, namely the QCD corrections to the electroweak p parameter which is given by 

p=l + 5p, (23) 

with 

, E z (0) S H ,-(0) 

5p= ^l~^- (24) 

Sy[/(0) and £z(0) are the transverse parts of the W and Z boson propagators defined 
through 

= 9j J^wizi (25) 

where are the corresponding polarization functions. The QCD corrections to Sp in 

the Standard Model (SM), where all quark masses except the top quark mass are set to 
zero, are known up to four loops [23-25,49-52]. In the following we consider a generic 
fourth generation of quarks which couples to the W and Z boson in the same way as the 
top and bottom quarks of the SM. For simplicity we denote the new doublet by (£', b'). 

£iy(0) can be obtained from C^_ l and via the relation 

W = I£f (26 » 

In the case of S^(0) there is no contribution form the vector part. Furthermore one has to 
consider additional contributions at three-loop order which are only present for currents 
where ipi = ip2 in Eq. (2). Thus we cast the Z-boson self energy in the form 

Ex(0) 



M% 



3G F m 2 t , ( - a /a s \2 (2) /a s \2 (2) \ 

T^Vl [ c -w* + {-) 6C -w> + IvJ J ' (27) 



where the quantities SC°^[ 2 l h and 5C , "i 2 s ing receive contributions from the double bubble 
and singlet Feynman diagrams involving two mass scales (see Fig. 9). The double bubble 
diagrams involving only one massive quark are contained in Ct.\ diag together with all 
other one-scale contributions. Note that the singlet-type contributions involving massless 
up, down, strange and charm quarks add up to zero. 

In Appendix A explicit analytical results are presented for diag , 5C°^[ 2 j b and SC^[^ . 
Let us mention that the treatment of 75 for the singlet diagrams can be found in Ref. [53]. 

Before presenting results we transform the MS quark masses m t > and my to the on-shell 
scheme [54, 55] since this is common practice in the context of analyzing electroweak 
precision data. For My = we reproduce the SM result [51,52]. The general result can 
be cast in the form 
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(a) (b) 

Figure 9: Additional contributions to the diagonal correlator with two mass scales: (a) 
double bubble and (b) singlet contribution. 

with X = My /Mf being the ratio of the on-shell quark masses. The two-loop result reads 
16 



A (D 



3 (1 + * 2 ) + ^ H (X) + 2X * (3 +J 4) H ,o(X) 



12 



(1 -X 2 



(29) 



which agrees with the findings of Ref. [56]. Using the asymptotically expanded master 
integrals we obtain a series expansion in X and 1 — X. The results including terms of 
order X 2 and (1 — X) 4 read 



A«(X) = * - 2» - ™*± + + ^ - U% + 

324 9 486 135 4 9 



4/4 32 
27 "Y 04 



486 135 
664C 3 



4 

27' 



+ 



972 6 4 



27 

4 11255 2 4 



1 13^_8Cj^ 2^ 
9 27 9 J 3 

, 4 ^2 4/ 4 1154/| 
9* k+ 27 7T 27 + ~ 9 



1943 



304/ 3 



A" 



32 



y a 4 + /^( - 



52 32tt 2 12* 



4G 



9 



tt z / 2 - 16tt z / 2 - — 



9 



+ 



4 
3 



57T 



A' 



52/ y 32/jy 



9 



9 / 9 9 

where lx = log X and = M 2 has been adopted 5 and 



3 

1651C 3 
27 



ni 



5 77T 2 

3 + 27 



(30) 



a 4 = Li 4 (l/2) « 0.51747906167389938633, 

S 2 = ^i=Im(Li 2 (e W3 )) « 0.26043413763216209896 



''Results for general /i can be found in rho.m [45]. 
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D 3 = 6C 3 -yC4-6[Im(Li 2 (e-/ 3 ))] 2 » 
Similarly, the expansion around My ~ M# leads to 



-3.0270094939876520198 . 



(31) 



A (2) (X) = (1-X) 2 



5933 16 



108 27 
116 , 8n t 



2807C 3 
72 



ni 



38 8tt s 



27 27 



9 



9 



465547 4tt 2 



4860 



+ 



+ 



135 135 



157781C 3 
2160 



ni 



473 2tt' 



810 135 



+ 0{{l-Xf). 



(32) 



In the file rho.m [45] analytical results up to order X 30 and (1 — X) 20 are provided. 

In the above equations n\ labels the number of massless quark flavours. In order to 
reproduce the SM result with massless bottom quark one has n\ = 4 and x = 0. The 
contribution of a fourth generation is obtained for ni = 6 which assumes a massless top 
quark at three-loop order. 

In the SM the higher order corrections to the p parameter are quite important. E.g., the 
three-loop corrections correspond to a shift in the top quark mass of about 2 GeV which 
is larger than the current experimental uncertainty. In the limit of equal quark masses Sp 
is proportional to the mass difference and thus the numerical impact of the higher order 
corrections is reduced as can be seen from Eq. (32). 



4 Summary 

In this paper moments of correlators formed by vector, axial-vector, scalar and pseudo- 
scalar currents are considered to three-loop order. The currents couple to quarks with 
different masses ni\ and m% which leads to two-scale vacuum integrals after expanding 
in the external momentum of the correlators. The moments can be expressed as a linear 
combination of 12 non-trivial three-loop master integrals which are discussed in details 
in Section 2. In particular, a Mathematica package TwoMassTadpoles .m is provided [45] 
which contains analytical results for all but six master integrals. For the latter high- 
precision numerical results are available. In a further Mathematica file coefhl2.m we 
also provide results for the first four moments for each correlator. 

As a by-product we compute three-loop corrections to the p parameter for a generic fourth 
generation of quarks with masses m t / and my for the up- and down-type quark. In addition 
to the calculation for the moments of the non-diagonal correlators one has to compute 
double-bubble and singlet-type diagrams where two different masses can be present in the 
fermion triangles. For the latter analytical results are provided; the complete results for 
the p parameter up to three loops is available in rho . m [45] . 
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A Appendix: Explicit results for C" ldi , c^C^'j 2 ^ and 



a,(2) 
sing 



In this appendix analytical results are provided for the building blocks of E^(0) in Eq. (27) 
using the MS definition of the quark masses. For convenience we choose for the renor- 
malization scale fi 2 = m 2 ; general results are available in electronic form from Ref. [45]. 
Defining 

-l,diag — — l,diag °-l,diag "r ^ J u -l,diag ' \°°) 



we have 



r»a>(2) 
°-l,diag 



K-? (i+ 4 + kI (i+i2)+3 ^ h ' w ) 

+ - ( - ^x 2 R {x) - 16x 2 H 0i0 (a;) - i- (l + x 2 ) (401 + 6vr 2 
e V 3 54 7 

\ 1 272 
-36Cs) ) + ^ 2 Ho(x) (-181 + 12tt 2 - 72C 3 ) - —x 2 R 0fi (x) 

+ ™x 2 RoMx) ~ ^ (1 + * 2 ) ( " 790 + (-5 + 40/ 2 ) tt 2 + yvr 4 - 40Z| 
II6OC3 - 960a 4 ) +n,| '"-^1 - "J^l + iil±f!l _ ^ Ho (,) 



2 H ,o,o(^) + ^(l + ^) (-l + 32C 3 )j. 



39 39 
+ Y * 2 H 0i0 (z) - yX 2 H , 0i0 (x) - 



The double bubble and singlet result is given by 

- <2) _ 4(l + x 2 ) i -1(1 + x 2 ) - lx 2 R (x) 



0U -l,db - n,3 ^ 
+ 



9e 3 e 2 
| (-1 + tt 2 ) (1 + a: 2 ) + |a: 2 H (x) + 16x 2 H , (x) 



e 

+ ~ (2 - (1 + tt 2 ) x 2 ) H (x) -Ul + 7x 2 ) H ,o(x) 
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2 (1 - 8x - 6x 2 - 8x 3 + x 4 ) H_ 1)0 , (x) 
3x 

256 2tt , , 2(l + 8x-6x 2 + 8x 3 + x 4 )Hioo(a;) 

-^(l + O (49 + 7T 2 - 38C 3 ) , 
^-Sng = &r (H 1)0j0 (x) + H_ lj0i0 (x)) -7(1 + x 2 ) ( 3 . (35) 

The three- loop result C°^[ 2 l iag depends on the number of massless flavours n\. Let us 
finally mention that the Matnematica package rho.m contains replacement rules which 
express the HPLs occurring in Eqs. (34) and (35) in terms of logarithms and dilogarithms. 
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